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ABSTRACT 



Recent observations with the Submillimeter Wave Astronomy Satellite 
indicate abundances of gaseous H2O and O2 in dense molecular clouds which 
are significantly lower than found in standard homogeneous chemistry models. 
We present here results for the thermal and chemical balance of inhomogeneous 
molecular clouds exposed to ultraviolet radiation in which the abundances of 
H2O and O2 are computed for various density distributions, radiation field 
strengths and geometries. It is found that an inhomogeneous density distribution 
lowers the column densities of H2O and O2 compared to the homogeneous case 
by more than an order of magnitude at the same Ay. O2 is particularly sensitive 
to the penetrating ultraviolet radiation, more so than H2O. The S 140 and p Oph 
clouds are studied as relevant test cases of star-forming and quiescent regions. 
The SWAS results of S 140 can be accommodated naturally in a clumpy model 
with mean density of 2 x 10^ cm~^ and enhancement /uv — 140 compared with 
the average interstellar radiation field, in agreement with observations of [C I] 
and (^^^CO of this cloud. Additional radiative transfer computations suggest 
that this diffuse H2O component is warm, ~ 60 — 90 K, and can account for 
the bulk of the lio — loi line emission observed by SWAS. The p Oph model 
yields consistent O2 abundances but too much H2O, even for [C]/[O]=0.94, if 
/uv < 10 respectively < 40 for a mean density of 10^ respectively 10^ cm^^. It 
is concluded that enhanced photodissociation in clumpy regions can explain the 
low H2O and O2 abundances and emissivities found in the large SWAS beam for 
extended molecular clouds, but that additional freeze-out of oxygen onto grains 
is needed in dense cold cores. 

subject headings: ISM: abundances — ISM: clouds — ISM: molecules 
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1. Introduction 

In contrast with the major carbon-bearing species (CO, C, C+), the abundances of the 
major oxygen-containing atoms and molecules in general molecular clouds are still poorly 
understood. The best-known species is CO, which reaches large abundances of ~ 10~^ when 
the shielding column of gas exceeds Ay ~1 for a typical interstellar cloud of ~ 10^ cm~^ 
illuminated by the average interstellar radiation field, locking up the bulk of the gas-phase 
carbon (e.g., Hollenbach & Tielens 1999, van Dishocck & Black 1988). Theoretical models 
indicate that at the edge atomic O is abundant, but that further into the cloud, more and 
more of the ultraviolet radiation is attenuated so that at Ay ~ 4 mag, CO ceases to be the 
dominant oxygen-bearing species and O2 reaches its peak abundance of > 10~^, with H2O 
present at the level of ~ 10~® (e-g-i Tielens & Hollenbach 1985, Lee et al. 1996). If true, 
then H2O and O2 could be major coolants in molecular clouds (Goldsmith & Langer 1978; 
Neufeld, Lepp & Melnick 1995). 

The launch of the Submillimeter Wave Astronomy Satellite (SWAS) (Melnick et al. 
2000) has made it possible to look for emission lines of H2O and O2 in cold molecular 
clouds. Although the Infrared Space Observatory (ISO) observed many II2O lines, it was 
sensitive primarily to warm, ~ 100 K molecular gas near young stellar objects where the 
H2O abundance can be enhanced due to ice evaporation and high-temperature (shock) 
reactions (e.g., van Dishoeck et al. 1999, Nisini et al. 2000). Measurements of, or upper 
limits on, the abundances of H2O and O2 provide crucial tests of many aspects of the 
oxygen chemistry in normal molecular clouds and of the cloud structure, and hence of our 
theoretical understanding of such systems. 

The initial SWAS results indicate that the abundances of H2O and O2 in dense 
molecular cloud cores are surprisingly low compared with the above model predictions, 
about 6 X 10^^° - 1 X 10"^ for H2O and < a few xlO"^ for O2 (Ashby et al. 2000, Snell et 
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al. 2000a, Goldsmith et al. 2000). In contrast, ISO observations of the [O I] 63 and 145 /xm 
hnes suggest substantial abundances of atomic O in quiescent gas (e.g., Baluteau et al. 1997, 
Caux et al. 1999). Bergin et al. (2000) propose a model in which freeze-out of oxygen on 
dust grains in the form of molecular ices is significant, with the remaining oxygen in atomic 
form in the gas, along with CO. While freeze-out and gas-grain interactions undoubtedly 
play a role in dense cold clouds, we present here an alternative explanation based on clumpy 
molecular clouds which may be more appropriate for the general lower density molecular 
cloud material contained in the large SWAS beam of 3.3' x 4.5'. The inhomogeneous nature 
of molecular clouds is well estabhshed from observations of extended [C I] and [C II] (e.g., 
Keene et al. 1995, Stutzki et al. 1988, Plume et al. 1999), and has been modeled by various 
groups (e.g., Meixner & Tielens 1993, Spaans 1996, Storzer et al. 1996). It is therefore 
natural to also explore the H2O and O2 abundances and emission in such models. 

2. Basic Model Description 

The results presented here were obtained by application of the numerical code of 
Spaans (1996), described further in Spaans & van Dishoeck (1997). The interested reader 
is referred to these papers for a description of the underlying algorithms. The code has 
been specifically designed to solve large chemical networks with a self-consistent treatment 
of the thermal balance for all heating and cooling processes known to be of importance 
in the interstellar medium (Spaans & Ehrenfreund 1999). The radiative transfer in the 
cooling lines is solved by means of a Monte Carlo approach with checks provided by an 
escape probability method for large line optical depths. The thermal balance includes 
photo-electric heating by dust grains, heating by cosmic rays, absorption of infrared photons 
by H2O molecules and their subsequent coUisional de-excitation, and gas-grain heating. 
The cooling includes the coupling between gas and dust grains (HoUenbach & McKee 1989), 
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atomic lines of all metals, molecular lines, and all major isotopes, as described in Spaans 
et al. (1994) and in Neufeld et al. (1995) for the regime of very large line optical depth. 
All level populations are computed in statistical equilibrium at a relative accuracy of no 
less than 10~^, an accuracy which is also imposed on the ambient line radiation field. The 
combined chemical and thermal balance is required to obey a convergence criterion between 
successive iterations of 0.5%. 

The adopted chemical network is based on the UMIST compilation (Millar, Farquhar 
& Willacy 1997), and is well suited for low temperature (<200 K) dense molecular 
clouds. A value of 0.3 is chosen for the branching ratio leading to H2O in the dissociative 
recombination of H3O+, with a rate coefficient of 3.3 x 10-^{T/300)-^-^ cm^ s^^ The 
branching ratio is consistent with the results of Vejby-Christensen et al. (1997) but higher 
than the value of 0.05 suggested by Williams et al. (1996). Observations and (homogeneous) 
models by Spaans et al. (1998) toward the translucent cloud HD 154368 favor a branching 
ratio of less than 0.3 (at 3a) so that the water abundance results presented here should be 
regarded as strict upper limits. The latest rate coefficients for important neutral-neutral 
reactions are adopted, in particular for C + O2 — CO + O and O + OH ^ O2 + H at 
2.48 X 10-12 (7^/300) i-54e-6i3/T ^^^^3 g-i ^^^d 1.77 x IQ-^^e'^^^/^ cm^ s-^ respectively. The 
reaction of O2 with neutral sulfur is not found to be a major sink for molecular oxygen. 
No gas- grain interactions are taken into account. The carbon and oxygen gas- phase 
fractions are set to 6c — 0.33 and 5o — 0.35 (C/O=0.45), but the case of Sq — 0.17 is also 
considered. These fractions are defined with respect to the solar values of [C] = 4.0 x lO"^ 
and [O] = 8.3 x lO""'. Depletions of other elements are as in Spaans & van Dishoeck (1997). 



-6- 



3. Model Results 

In order to investigate the dependence of O2 and H2O on density, radiation field 
and geometry, spherical and slab-like model clouds were computed for homogeneous and 
inhomogeneous density distributions (Spaans 1996). The clumpy models are defined by 
two parameters: the volume filling factor /, which fixes the fraction of the total volume 
that is occupied by the clumps in a two-phase density medium, and a characteristic clump 
size ^c) which fixes the extinction through an individual clump. The ratio between the 
high (h) density and low (1) density phase is chosen to be r = 20. It then follows that 
the mean total hydrogen density n obeys n — fn\i -\- (1 — f)n\. We model two types of 
regions: a star-forming cloud where the radiation field is significantly enhanced due to a 
nearby O or B star (e.g., S 140) and a more quiescent region such as p Oph. For the S 140 
extended molecular cloud, Spaans & van Dishoeck (1997) have constrained / <50% and 
4 > 0.2 pc fi-om observations of [C I] and ^^CO (Plume et al. 1994), usmg an average 
density n ^ 10^ cm^^. Recent ISO data of [C II] and [O I] confirm densities of order 10^ 
cm~^ in the extended cloud (Li et al. 2001). We adopt here / — 30% and £c — 0.4 pc with 
n = 2 X 10^ cm~^. For models with a larger mean density the clump size Ic is decreased by 
the same factor with which the mean density is increased. The enhancement of the radiation 
field is taken to be /uv=140 with respect to the Draine (1978) field, consistent with the 
enhancement at the edge of the cloud near the BOV star HD 211880. This enhancement 
factor decreases to /uv ~ 30 — 50 further into the extended cloud due to geometric dilution. 
The average interstellar radiation field incident from other directions has been included as 
well. The p Oph cloud has been studied by several groups (see Liseau et al. 1999) and is 
modeled with / = 20% and r = 30, more clumpy than S 140 but with an identical value for 
Ic- We adopt /uv 10 and n — lQ^ cm~^. It thus provides an interesting comparison. 

Figures 1-2 present the basic results of this work. For the inhomogeneous models the 
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shell/slice averaged abundances are shown for the sphere/slab as functions of Ay. It was 
found that the dispersion around the mean values by taking different cuts is not larger than 
about 50%, so that the systematics discussed below are robust. The following trends can 
be identified. First, for the same Ijjy and n, a strong decrease in the abundances of O2 and 
H2O is found in the inhomogeneous models, although less so for H2O. The reason is that 
the destruction of O2 remains dominated by photodissociation to large extinctions and that 
the clumpy structure allows ultraviolet photons to penetrate to larger depths. The effect 
spans a factor of 30 or more. Second, the central abundances of H2O and O2 are decreased 
by a factor of 2-4 as the geometry is changed from a plane-parallel slab to a sphere. The 
reason is simply that the shielding column is larger for some rays emanating from the 
center of a slab compared to a sphere. A clumpy density distribution is quite effective in 
redistributing the directed ultraviolet radiation field from a (point-like) source into a more 
isotropic distribution at a depth of a few magnitudes of visual extinction (Spaans 1996). 
Third, a systematic trend is observed with varying dissociation parameter U — luy/n. 
When this number is small, < 5 x 10"^, i.e., at low /uv and/or high n, H2O and O2 rise 
sharply at the edge of the cloud, and H2O quickly reaches the regime where its removal is 
dominated by chemistry rather than photodissociation. The total observed column of H2O 
and O2 varies greatly with total extinction along the line of sight and it would be difficult 
to distinguish between homogeneous large U and inhomogeneous small U situations based 
on H2O and O2 alone. Of course, large U regions are bright in lines such as [C II] 158 //.m, 
[O I] 63 and 145 /im and high- J CO, whereas small U regions tend to have more prominent 
[C I] 609 /im emission, in particular when they are clumpy. 



4. Discussion and Compcirison with Observations 
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4.1. H2O and O2 abundances 

The observational results represent abundances integrated over depth. Thus, mass- 
weighted abundances for the spherical models have been computed for comparison with 
observations. Specifically, the mean abundance X{r) enclosed within some fractional radius 
r is given by integration of the abundance A{r) through X{r)—jQ A{r')r''^dr' / jQr''^dr', 
where is at the center and 1 at the edge of the model cloud. Note that only 1/8 of the 
mass is contained inside the half radius of the spherical constant density model cloud, while 
this would be 1/2 for a 1/r^ density profile. 

For the clumpy, spherical model of Figure 1 (top and middle) the above approach leads 
to a predicted mean H2O abundance ^(1) = 3 x 10~® and an O2 abundance of 4 x 10~^° 
for a typical total depth of Ay — 10 mag. Snell et al. (2000a) detected H2O emission 
from the S 140 dense star-forming core and PDR region with SWAS, but not from the 
surrounding extended molecular cloud. Their inferred abundances of ~ 10~^ (see also the 
detailed analysis by Ashby et al. 2000) for the core and < 10~^ for the extended cloud are 
only a factor of three lower than our model results and consistent within the factor of three 
or more uncertainty in their assumed density (the inferred H2O abundance scales with the 
inverse of density). Figure Ic (bottom) shows the effect of increasing the density to 5 x 10"^ 
cm~^ for the S 140 environment. The smaller value of U causes photodissociation to play a 
lesser role in the removal of O2 and H2O, and hence dumpiness is of less importance. The 
model (mass weighted) abundances are 3 x 10~^ and 4 x 10~^ for H2O and O2, respectively. 

The model for p Oph appears less successful in reproducing the observed abundances 
of O2 and H2O. The mass weighted H2O abundance for Figure 2c (bottom) is 3 x 10~^, 
a factor of ~100 above the measured value. Panels 2a and b (top and middle), which are 
representative of larger values of U and roughly consistent with the results of Liseau et al. 
(1999) who found nn = 3 x 10^ — 3 x 10^ cm^^, yield mass-weighted H2O abundances of 
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5 X 10~^ and 2 x 10~^, respectively. Prom the discussion in §3, it follows that the S 140 case 
yields better agreement with the observations mainly because of its larger illumination (and 
U parameter). The clumpy models enhance the importance of photodissociation of H2O 
and O2 deeper into the cloud, and therefore work best for a strong incident radiation field. 

In all inhomogeneous models with n ~ 10^ cm~^, the mass weighted abundance of O2 
is below 10~^, both for S 140 and p Oph, consistent with the overall SWAS upper hmits as 
well as the specific upper bounds for S 140, < 7 x 10~^, and p Oph, < 3 x 10"^ (Goldsmith 
et al. 2000). Only the quiescent model of Figure 2c (bottom) with high density and low 
IjjY, yields a mass weighted O2 abundance of 10~^. This is a factor of a few above the 
observational hmit (Goldsmith et al. 2000), but this is easily alleviated by a modest increase 
in the [C]/[0] ratio (Bergin et al. 2000). To test the latter possibility, we have run a model 
with an oxygen depletion of 0.17 and [C]/[O]=0.94. A decrease in the mass weighted H2O 
(O2) abundance of a factor of 5 (20) is found, bringing the model in agreement with the 
observations for S 140 and factors of 20, 4 and 2 too high for H2O toward p Oph in models 
2a, b and c, respectively, although a smaller value for the HaO"^ branching ratio of 0.05 
instead of 0.3 could further reduce these discrepancies. 

4.2. H2O and O2 emission 

In order to verify whether the low abundances of H2O and O2 found for the clumpy 
model clouds can also reproduce the bulk of the observed emission, radiative transfer 
computations have been performed using a Monte Carlo method (Spaans 1996). It was 
found that two effects dominate the emissivity of the H2O lio — loi line: temperature and 
geometry. The temperature is very important because it influences both the ortho-para ratio 
of the main coUision partner H2, 0.1 at 30 K and 1.0 at 80 K (Neufeld & Sternberg 1999), 
as well as the coUisional excitation rate. The H2 ortho-para ratio is a crucial parameter 
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to compute rigorously from the thermal balance and FUV pumping since the 0-H2 J—1 
collision rates with H2O are more than an order of magnitude larger than the P-H2 J=0 
rates (Phillips, Maluendes & Green 1996). The geometry is important because the lio — loi 
hne is optically thick but effectively thin. This implies that line photons emitted by the 
warm edges of a spherical clump can interact with water molecules in the colder interior, 
and so constitute an excitation term that raises the level population in the lio state. Note 
in this that the curvature of the sphere causes the warm gas to cover 47r steridians, unlike 
the slab case. It is found that a spherical geometry, for the S 140 temperature profile 
given below, raises the lio — loi line strength with a factor of 2.4 compared to that from 
the plane-parallel slab. Furthermore, for identical physical conditions, a plane-parallel 
slab model yields a temperature profile that is systematically 20-30% cooler than for the 
corresponding spherical case, mainly because the latter possesses more lower extinction 
lines of sight. 

The thermal balance yields a temperature of 150 (80) K at the edge of S 140 (p Oph) 
and 20 (10) K at its center. The mass weighted H2O temperatures are ~ 80 K and ~ 30 K 
for S 140 and p Oph, respectively. For S 140, these values refer to the model of Fig. la/6, 
and for p Oph to the model of Fig. 2b. The S 140 value is significantly larger than the typical 
temperature in dense cores, ~ 30 K. The resulting integrated hne intensities convolved with 
the SWAS beam are 1.6 K km s~^ and 3.8 K km s~^ for S 140 and p Oph, respectively. 
The clumpy S 140 model can thus reproduce the observed value (1.7 K km s""^) quite 
nicely, but the clumpy p Oph H2O emission, just as the abundance, is too large compared 
to the observed value of 0.8 K km s~^. Use of a plane-parallel geometry would reduce the 
discrepancy for p Oph significantly. In both instances, the bulk of the signal in the SWAS 
beam is associated with the PDR. In fact, the H2O emission arises mainly from the warm, 
60-90 K, edges of irradiated clumps. Therefore, a strong correlation should exist between 
the SWAS H2O emission and high— J (e.g. 5-4, 6-5) CO emission. Such a correlation is in 
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fact observed toward M 17SW (Snell et al. 2000b). It is the elevated temperatures at the 
surfaces of the clumps that makes the results so sensitive to the ortho-para ratio. For the 
case of S 140, the extended molecular cloud gives an emission that is a factor of 30 weaker 
at a position 4' north-east of the PDR, due to the lower ambient temperature, consistent 
with observations. The O2 Nj = 83 — 12 line is truly optically thin and yields, due to its 
simpler excitation, more robust results. We find 1.8 x 10^^ K km s~^ and 2.3 x 10^^ K km 
s~^ for the S 140 and p Oph case, respectively. 

ISO H2O observations suggest that much of II2O toward star-forming cores such 
as S 140 is associated with a small region with enhanced H2O surrounding the young 
stellar object (s) (e.g., Wright et al. 1997, Ceccarelli et al. 1999). Since this emission is 
subsequently beam-diluted in the large SWAS beam, its contribution could be small and 
still be consistent with the fact that the lower density gas in the SWAS beam is capable 
of reproducing the bulk of the detected water emission by SWAS. We have verified this 
by adding a clumpy star-forming core with n = 10^ cm~^ and a size of ~ 0.5 pc (Zhou 
et al. 1994), and the same values for / and r, to the S 140 extended molecular cloud as 
in Spaans & van Dishoeck (1997). Even though the mean water abundance goes up by 
a factor of 30, the lower kinetic temperature and smaller spatial extent yield a relative 
contribution of ~ 20% after convolution with the SWAS beam. Furthermore, the Zhou et 
al. CS observations can be well reproduced by this clumpy core that includes densities of 
up to ~ 10^ cm^^ and temperatures that are still ~ 30 — 40 K. Q 

Bergin et al. (2000) have argued against the importance of photodissociation because 

^At even higher densities (~ 10^ cm~'^) and smaller scales (~ 200 — 400 AU), the H2O 
abundance becomes even larger and the region emits thermally at about 100 K. Strong 
beam-dilution limits this contribution to the SWAS signal to less than 2% at the distance of 
~ 910 pc for S 140. 
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it would affect the abundances of other species such as NH3 as well. The S 140 model, 
including the star-forming core as above, was checked for consistency with the observational 
NH3 values of 2 — 9 X 10^"^ cm~^ (Ungerechts, Winnewisser & Walmsley 1986). Our model 
value of ~ 4 X 10^^ cm"^ is not in conflict with observations. In addition, even a moderate 
enhancement in NH3 formation due to grain-surface chemistry is likely to resolve any mild 
discrepancy. For the extended molecular cloud outside the core there are no observational 
constraints for ammonia. Future observational and theoretical studies should investigate 
other species whose abundances are sensitive to photodissociation, such as the CN/HCN 
and C2H/HCO"'" ratios (e.g., Jansen et al. 1995). This work has not addressed the effects 
of time-dependent chemistry and freeze-out on the O2 and H2O abundances (Bergin et al. 
2000, Viti 2000, private communication). Our results suggest that dumpiness is a viable 
alternative explanation for extended lower-density molecular clouds in star-forming regions 
like S 140 with high radiation field, but they also strengthen the case for time- dependent 
chemistry and freeze-out for a more quiescent dense region like p Oph. The freeze-out 
models may in fact benefit from adoption of a clumpy structure since the central dust 
temperature is enhanced by the higher internal ultraviolet field, thus providing a mechanism 
to moderate the very rapid freeze-out. Future observatories such as ODIN and the HIFI 
instrument on the Far Infrared and Submillimeter Telescope (FIRST) will have increased 
sensitivity, and, in the case of FIRST, a much smaller beam. Therefore, more definite tests 
of the various classes of models will be forthcoming in this decennium. 

We are grateful to Ted Bergin, Gary Melnick, Matt Ashby and Serena Viti for 
stimulating discussions regarding the SWAS results, to Steven Doty for assistance with the 
chemistry, and to the referee, Ted Bergin, for his constructive comments that helped to 
improve this paper. 
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Fig. 1. — Models appropriate of the S 140 extended molecular cloud, la (top): Model cloud 
with n — 2 X 10^ cm~^ and /uv = 140. Sohd lines are clumpy spherical model, dashed 
lines homogeneous spherical model; lb (middle): Model cloud with n = 2 x 10^ cm~^ and 
7uv — 140. Solid lines are clumpy spherical model, dashed lines clumpy slab model; Ic 
(bottom): Model cloud with n — b x 10^ cm~^ and /uv = 140. Solid lines are clumpy 
spherical model, dashed hues homogeneous slab model. 



Fig. 2. — Models appropriate for the p Oph cloud. 2a (top): cloud with n — 10^ cm~^ and 
7uv = 10. Solid lines are clumpy spherical model, dashed lines homogeneous slab model; 2b 
(middle): Model cloud with n = 10^ cm~^ and Juv = 40. Solid lines are clumpy spherical 
model, dashed hues homogeneous spherical model; 2c (bottom): Model cloud with n — 
cm"^ and /uv — 10. Solid lines are clumpy spherical model, dashed lines homogeneous 
spherical model. 



